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ABSTRACT 

Eight dwarf starburst galaxies have been observed with the Caltcch SubmiUimeter Observa- 
tory (CSO) telescope in the CO J= 3-2 transition. The galaxies observed are He 2-10, NGC 
5253, NGC 1569, NGC 3077, Haro 2, Haro 3, II Zw 40 and Mrk 86; all but the last two are 
detected. The central regions of He 2-10 and NGC 5253 were mapped and a C0(2-l) spectrum 
of NGC 5253 was obtained. The error weighted mean CO(3-2)/CO(l-0) ratio of the detected 
galaxies is 0.60±0.06, which is virtually identical to what is found for starbursts in the nuclei 
of nearby spirals, and suggests that the molecular gas is optically thick, warm (Tfe > 20 K), 
and moderately dense (n^r^ ^ 10"^"^ cm~^). The CO(3-2)/CO(l-0) ratio peaks at or close to the 
starburst in all cases. CO emission does not appear to be optically thin in these dwarfs, despite 
the low metallicity and intense radiation fields, which is probably because in order for CO to exist 
in detectable amounts it must be self-shielding and hence optically thick. Physical properties of 
the molecular clouds in these dwarf starbursts appear to be essentially the same as nearby spiral 
nuclei, with the possible exception that CO is more confined to the cloud cores. 

Subject headings: galaxies: dwarf — galaxies:ISM — galaxies: nuclei — galaxies: starburst — radio 
lines: galaxies 



1. Introduction 

Dwarf galaxies differ in many characteristics 
from large star-forming spirals. Typically they are 
much smaller (<3 kpc), have lower rotation veloc- 
ities (<30 km/s), generally warmer ISMs (T^; > 
30K), solid body rotation, and lack density waves, 
presenting less complicated star forming environ- 
ments than spiral galaxies (eg. Hodge 1971; Gal- 
lagher and Hunter 1984; Thronson and Telesco 
1986). Also, dwarf galaxies are compositionally 
different, appearing to have relatively high HI and 
dark matter fractions (Thuan and Martin 1981), 
relatively little dust, and lower metallicities. The 
molecular content in dwarf galaxies is particularly 



difficult to study, since CO (1-0) is usually quite 
weak due to the low metallicities, and therefore 
probably not a good tracer of H2 (Israel et al. 
1986; Verier and Hodge 1995; Wilson 1995). 

The relatively weak CO(l-O) emission in dwarf 
galaxies leads us to search for a stronger line 
with which to study molecular gas. The J=3- 
2 transition of CO may be more easily detected 
in dwarf starbursts since it is relatively easily ex- 
cited in warm gas (E/k = 33.2 K). Due to its 
higher characteristic temperature and critical den- 
sity {ricr ~ 2 X 10^ cm~^), the CO(3-2) transition 
may be more sensitive to warm, dense gas directly 
involved in the starburst. It also has a higher op- 



tical depth than the CO (1-0), which can make it 
easier to detect if CO (1-0) is optically thin. 

A sample of nearby dwarf starbursts have been 
observed in the CO(3-2) transition with the Cal- 
tech Submillimeter Observatory (CSO). Among 
the questions addressed are: do molecular clouds 
in dwarf galaxies have different physical proper- 
ties than their higher metallicity counterparts? Is 
the CO(l-O) line weak in dwarf galaxies because 
it is optically thin? Is the CO conversion factor 
in dwarf starbursts different from large spirals be- 
cause of these physical conditions? 

1.1. The Galaxy Sample 

A sample of eight nearby {^20 Mpc) dwarf star- 
burst galaxies was observed. The galaxies are He 
2-10, NGC 5253, NGC 1569, NGC 3077, Haro 
2, Haro 3, II Zw 40 and Mrk 86. The sample 
is heterogeneous and illustrative, but not com- 
plete or unbiased. Galaxies were selected based 
on the following criteria: they have been previ- 
ously observed and detected in either the CO(l-O) 
or C0(2-l) transition (Wikhnd and Henkel 1989; 
Becker, Schilke and Henkel 1989; Sage et al. 1992; 
Baas, Israel and Koornneef 1994; Greve et al. 
1996), have IRAS 100/im fluxes greater than 5 
Jy (Thronson and Telesco 1986; Melisse and Is- 
rael 1994), are fainter than Mb > —18.5 and 
have undergone a recent burst of intense star- 
formation. All the galaxies except NGC 3077 and 
Mrk 86 have Wolf-Rayet (WR) emission features 
(eg. Conti 1991; Steel et al. 1996; Gonzalez Del- 
gado et al. 1997; Schaerer, Contini and Pindao 
1999). Some of these galaxies contain super star- 
clusters (SSCs) (He 2-10, NGC 5253, NGC 1569 
and Haro 3) (Conti and Vacca 1994; Meurer et al. 
1995; Arp and Bandage 1985; Steel et al. 1996). 
A number of these dwarfs show signs of interac- 
tion (He 2-10, NGC 5253, NGC 1569, NGC 3077 
and II Zw 40). The properties of the galaxies in 
the sample are shown in Table 1. 

2. Observations 

We observed the ^^CO{3-2) line (345.796 GHz) 
towards the galaxies He 2-10, NGC 5253 and 11 
Zw 40, on 1997 February 24 & February 25, and 
the remaining five galaxies on 1999 January 14 
& January 15 using the 10.4m Caltech Submil- 
limeter Observatory (CSO). The beamsize of the 



CSO at 345 GHz is 22". Facihty SIS receivers 
were used together with a 1024 channel 500 MHz 
AOS spectrometer. System temperatures ranged 
from 600 K - 1100 K for the 1997 data and 700 
K - 900 K in 1999 (T225 ^ 0.08 - 0.15). Point- 
ing was checked using IRC-M0216, a Ori, (1997 
observations) Saturn, and Mars (1999). Absolute 
pointing uncertainty was <6 , and repeatable to 
< 3 — 4 . Reported temperatures are main beam 
temperatures, which is the brightness temperature 
a source would have if it uniformly filled the main 
beam and was zero elsewhere. The main beam 
efficiencies used to convert the antenna tempera- 
ture to the main beam temperature were deter- 
mined from observations of Mars and found to be 
Vmb = 0.6 (1997) and 0.62 (1999). Since no pub- 
lished C0(2-l) spectrum exists for NGC 5253, G. 
Serabyn kindly obtained one for us on 1997 May 
28. At 230 GHz, the CSO beamsize is 30", and 
Vmb = 0.65 as determined from observations of 
Mars. System temperatures for this spectrum are 
~400 K. Second order or lower polynomials were 
removed from the each scan to give a flat base- 
lines. The scans were then averaged, using l/T^yg 
weights. The spectra were Hanning smoothed to 
a resolution of 8.0 km s~^ (1997) and 6.6 km s~^ 
(1999). Reduction of the spectra were done in 
CLASS. For He 2-10, the maps were made and 
analyzed using the NRAO AIPS package. 

2.1. The CO(3-2) Spectra 

CO(3-2) spectra of the central positions of each 
of the galaxies are displayed in Figure 1. We de- 
tected CO(3-2) in all the galaxies except Mrk 86 
and II Zw 40. Line intensities, imcertainties, and 
gaussian fits are tabulated in Table 2. For some of 
the spectra gaussians may not represent the true 
line profile but given the weakness of the signal, 
fitting anything more complicated is unwarranted. 
The details of each galaxy are discussed in §6. 

In general CO (3-2) is very weak (Tmb ^ 0.05 
K) over our projected beamsizes of 0.25 - 2 kpc. 
This is not surprising; dwarf galaxies are known 
to be weak in the other low J CO lines (Young, 
Gallagher and Hunter 1984; Israel 1986; Tacconi 
and Young 1987; Arnault et al. 1988; Sage et 
al. 1992). However, CO (3-2) is strong relative 
to CO(l-O) and C0(2-l). Only one of the six de- 
tected galaxies has CO(3-2)/CO(l-0) ratio as low 
as 0.4, which is typical of the Galactic disk. Com- 
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Fig. 1. The CO (3-2) spectra of the eight galax- 
ies. The vertical axis is main-beam temperature 
in units of K. The horizontal axis is velocity (LSR) 
in units of km s"^. The line represents the least- 
squares fitted gaussian. 



paring the CO(3-2)/CO(l-0) line ratio with the 
CO(2-1)/CO(1-0) line ratio, we find that galaxies 
with high 2-1/1-0 ratios tend to have high 3-2/1-0 
ratios. There is one notable exception, II Zw 40. 
II Zw 40 has a high CO(2-1)/CO(1-0) ratio (Sage 
et al. 1992), but its CO(3-2)/CO(l-0) ratio is 
<0.43 {la). 

3. Line Ratios: Physical Properties of the 
Molecular Gas 

CO line ratios can be used to investigate molec- 
ular gas properties in these starbursts. C0(3- 
2)/CO(l-0) ratios are more sensitive to gas tem- 
perature than CO(2-1)/CO(1-0) ratios because of 
the larger separation in J. Thanks to the fortu- 
itous matching of the CO(l-O) beam at the IRAM 
30m (~2l") and the CO(3-2) beam at the CSO 
(~22' ), CO(3-2)/CO(l-0) line ratios can be esti- 
mated without assumptions about source struc- 
ture. Only the two southern galaxies, He 2-10 
and NGC 5253, do not have IRAM CO(l-O) spec- 
tra. For these two galaxies, the source has been 
mapped at CO (3-2) or with an interferometer at 
C0(2-l) or both (Meier and Turner 1998). 

Along with the CO(3-2)/CO(l-0) line ratios. 
Table 3 includes CO(2-1)/CO(1-0) line ratios. 
The C0(2-l) observations do not have match- 
ing beams. Therefore Table 3 indicates what ra- 
tio would be predicted in the point-source limit 
(source much smaller than either beam) and in 
the uniform filling limit (source much larger than 
both beams). These line ratios are calculated us- 
ing (eg. Rohlfs and Wilson 1996, page 362): 

where fi^, "^^Hb and are the source solid an- 

gle and the solid angle of the beams at the re- 
spective transitions, and I21 and /lo are their ob- 
served line intensities. For sources which have 
been mapped a deconvolved ratio based on the 
measured Clg is also included. In all cases, with 
the possible exception of the CO(2-1)/CO(1-0) in 
NGC 1569, the uniform filling ratios are unrealis- 
tically high, implying the sources are small. This 
is consistent with the low filling factor found for 
these galaxies (§3.1) 

The CO(3-2)/CO(l-0) ratios in the dwarf star- 
bursts range from 0.37 - 1.1. The error weighted 
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mean CO(3-2)/CO(l-0) line ratio for the sample is 
0.6 ±0.06. (We have estimated the uncertainty in 
the line ratio as the noise error in each spectrum 
added in quadrature with 20% absolute calibration 
uncertainties.) This error weighted ratio is iden- 
tical, within the uncertainties, to what is found 
for a sample of non-dwarf starburst nuclei (0.64; 
Devereux et al. 1994) and for nearby luminous 
IR galaxies (~0.7; Mauersberger et al. 1999). The 
value is higher than the value obtained for Galac- 
tic GMCs, 0.4, and closer to the value obtained 
for star- forming cores in the Galaxy (0.6; Sanders 
et al. 1993). 

3.1. LTE Modeling: Optically Thick CO 
Emission 

Comparing CO (3-2) to the other CO transitions 
(CO(3-2)/CO(l-0) line ratio is used unless other- 
wise stated) can, in principle, be used to constrain 
gas excitation temperatures and optical depths. 
The ratio of the CO(3-2) intensity to the CO(l-O) 
intensity, under the LTE assumption is: 

J32 _ / T32dv 
ho ~ jTiodv 

^ hl{Jsi{Tex) - ■/32(rcmb))(l - e-'^) 
fw{Jw{Tex) - JwiTcmbW - e"^!") 

where J,ATex) {hv / k) / {exp{hv / kT^-jc] - 1), 

with T and / being the optical depth and filling 
factor of each transition, respectively (eg. Hurt et 
al. 1993). To increase the signal-to- noise ratio 
(SNR), we have used the ratio of integrated in- 
tensities instead of the peak main-beam tempera- 
tures. This requires that the CO(3-2) and CO(l-O) 
line profiles be similar. For most of the galaxies 
this is a good approximation (§6). For the LTE 
approximation, ratios greater than unity indicate 
warm, optically thin gas, while ratios less than 
unity indicate optically thick molecular gas. 

In practice, non-LTE effects such as different 
source sizes for the CO (3-2) and CO (1-0) tran- 
sitions or different Tex for each transition due 
to temperature gradients and differential optical 
depth effects in externally heated clouds may af- 
fect the ratio. This dataset cannot address such 
details. A discussion of these details would require 
higher resolution or observations of rare CO iso- 
topomers. The ^^CO isotopomers, whose transis- 
tions will be less optically thick, can be helpful in 



better constraining the properties of these galax- 
ies, but unfortunately these galaxies are too weak 
in ^^CO to detect with current telescopes (except 
He 2-10). 

Previous studies indicate that variations in Tex 
become important in localized regions of high col- 
umn density, where optical depths of order unity 
are reached over very small physical distances (eg. 
Turner, Hurt and Hudson 1993; Meier, Turner 
and Hurt 2000). When averaged over large single- 
dish beams, the potential error associated with 
such effects appear to be minor. (For example, 
in the case of the nucleus of the nearby metal-rich 
starburst, IC 342, a comparison of the physical 
conditions derived using single-dish ^^CO observa- 
tions (Eckart et al. 1990) are similar to those ob- 
tained with higher resolution, even with the pres- 
ence of strong local temperature gradients (Meier, 
Turner and Hurt 2000). Similar results have been 
found from studies of Galactic star-forming re- 
gions (eg. Wilson, Howe and Balogh 1999). 

Nearly all of the galaxies we observe have ratios 
$1.0 which indicate optically thick emission. NGC 
3077 (1.1) and the one velocity component of Haro 
3 (1.3) have ratios that are larger than one, but 
only marginally so. In the case of NGC 1569 and 
NGC 3077, some portion of the high line ratios are 
due to a slightly larger linewidth of the CO (3-2) 
line relative to the CO (1-0) lines, indicating that 
there maybe some regions off line center that have 
CO(3-2)/CO(l-0) ratios greater than 1.0. But in 
both cases the SNR is too low to say with cer- 
tainty. It appears that optically thick gas domi- 
nates the CO emission in these dwarf galaxies. 

We estimate the gas excitation temperatures of 
the six detected galaxies, assuming the molecular 
gas is optically thick. Derived Tex range from ~ 6 
K (Haro 2) to >50 K (NGC 1569, NGC 3077, Haro 
3). For the high ratio galaxies a specific excita- 
tion temperature cannot be obtained because the 
CO(3-2)/CO(l-0) ratio ceases to differentiate tem- 
peratures well for Tft:<30 K. As a result of beam 
dilution, Tmb/Tcx gives an estimate of the areal 
filling factor for each of the galaxies. Rather low 
filling factors are found, with values ranging from 
fa ~0.02 to <0.001 (Table 4). 
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Fig. 3.— The CO(3-2) intensity in K km s"^ is 
plotted versus metallicity. Bold open squares are 

the CO (3-2) detections, and bold crosses arc up- 
per limits. For reference, the data is plotted over 
the dwarf galaxy CO(l-O) spectra compiled from 
the literature by of Taylor ct al. (1998). Diamonds 
represent detections and the arrows are upper lim- 
its. 



Fig. 2. — The LVG models for the six detected 

galaxies. The solid line is for the C0(2-1)/C0(l- 
0) line ratio and the dashed line is for the C0(3- 
2)/CO(l-0) line ratio. For each case, the bold line 
represents the measured value and the thin lines 
represent the ±la range. The best fit solutions 
are the regions where both line ratios overlap. All 
models displayed here use a velocity gradient of 1 
km s~^pc~^. The [C0/-ff2] abundances are 8 x 
lO-'"^ (1 Zq) for He 2-10 and NGC 3077, 2.7 x 
10-5 (1/3 Zq) for Haro 2 and Haro 3, and 1.1 x 
10-5 (1/7 Zq) for NGC 1569 and NGC 5253. For 
He 2-10, the i2cO(l-0)/i3co(l-0) line ratio was 
also modeled. The grey curve represents the lower 
limit of Baas et al. 1994 over the same beamsize 
as the CO(3-2) line ratios (and consistent with the 
value obtained by Kobulnicky et al. (1995) over a 
larger beamsize). The abundance per unit velocity 
gradient used in this model is 5 x lO-'' km s-^pc-^ . 



o 
o 
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Beam Size (Kpc) 

Fig. 4.— The CO(3-2)/CO(l-0) line ratio, as a 
function of the projected physical size of the beam 
over which the ratio has been measured, for the 
galaxies detected from our sample. A weak corre- 
lation is seen. 
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3.2. LVG Modeling: Molecular Cloud 

Densities 

Additional constraints on the physical condi- 
tions of the molecular gas can be obtained us- 
ing a Large Velocity Gradient (LVG) model (Gol- 
dreich and Kwan 1974; De Jong, Chu and Dal- 
garno 1975). Due to the relatively large Doppler 
shifts within the cloud, the emergent intensity can 
be related to its local sources of excitation, 
and Nh^. As a result, LVG models can be used 
to get an indication of the gas kinetic temper- 
ature and density, given observations of several 
line ratios or antenna temperatures. While fill- 
ing factors do not strongly effect the line ratios 
for our matched beam observations, they do sen- 
sitively effect the brightness temperatures. The 
molecular clouds in these galaxies arc unresolved. 
Therefore, only line ratios provide useful con- 
straints. Detailed LVG modeling is not possible 
with the two line ratios, CO(3-2)/CO(l-0) and 
CO(2-1)/CO(1-0), since there are three unknowns, 
the gas density, ur^ , the kinetic temperature, Tk 
and Xco / dv / dr. As mentioned earlier, tempera- 
ture gradients can also contribute; however, for a 
reasonable assumption of the abundance per ve- 
locity gradient it is possible to provide rough con- 
straints on the physical parameters. Given our 
limited knowledge of the physical conditions in 
dwarf starbursts, even rough constraints are im- 
portant. 

The LVG models were run for gas densities 
ranging from njy^ = 10^ — 10® cm~^ and kinetic 
temperatures Tk 1 — 100 K, for a CO abundance 
per velocity gradient, Xco/dv/dr, ranging be- 
tween 10~® and 10""^. Given that we only have two 
lines to constrain three variables, model solutions 
are displayed for only one choice of Xco/dv/dr. 
Figure 2 displays the best fit solutions for Tk 
and with an Xco/dv/dr estimated in the 
following manner. The CO abundance is as- 
sumed to be the solar CO abundance scaled rel- 
ative to the galaxies' metallicity and a constant 
velocity gradient of 1 km s~^pc~^, which should 
roughly account for metallicity effects (Verter and 
Hodge 1995; Wilson 1995). The solutions ob- 
tained are relatively insensitive to Xco/dv/dr as 
long as it is within a factor of 4-5 times the dis- 
played value. For different assumed values of 
Xco/dv/dr, the best fit solutions shift shglitly to- 
wards higher (lower) densities, for lower (higher) 



values of Xco/dv/dr, due to the rise of the effec- 
tive critical density caused by decreased radiative 
trapping. 

Sage et al. (1992) argue that C0(2-l) /CO(l-O) 
line ratios obtained by approximating the source 
size as a point source are closer to reality because 
of the small CO source sizes expected for dwarf 
galaxies. The low derived CO(3-2) filling factors 
appear to confirm this. Therefore, the "point- 
source" ratio approximations are used for galax- 
ies with no mapping information (Table 3). The 
details of the LVG solutions for each galaxy are 
ad(lrc;ssed in §6, but the basic results can be un- 
derstood in simple terms. In order to excite a 
bright CO (3-2) line the gas must either be warm 
or have a relatively high density. At high densities 
the lower J transitions of CO will be thermalized 
and the CO(2-1)/CO(1-0) and CO(3-2)/CO(l-0) 
values will just reflect the Rayleigh-Jeans temper- 
ature corrections. At low densities, CO(3-2) weak- 
ens faster than C0(2-l) because of its higher criti- 
cal density. Therefore, galaxies with similar C0(2- 

1) /CO(1-0) and CO(3-2)/CO(l-0) values indicate 
thermalized, high density gas (the larger the ra- 
tios the higher the temperature) , while in galaxies 
with CO(3-2)/CO(l-0) < CO(2-1)/CO(1-0) lower 
density gas is indicated. 

The galaxies with low to moderate C0(3- 

2) /CO(l-0) and CO(2-1)/CO(1-0) values (He 2- 
10, Haro 2 and Mrk 86) are best fit by molecular 
gas that is dense (n//^ > 10** cm~^) and cold 
(Tk > 5 - 10 K). The galaxies with high C0(3- 
2)/CO(l-0) and CO(2-1)/CO(1-0) values (NGC 
1569, NGC 3077 and NGC 5253) have a large 
range of parameter space that fit the observed 
values. There is a warm, low density solution 
{Tk > 40 K; hh^ - lO^ '^ cm.-'^; NGC 5253) or a 
warm, high density solution (n/f^ > lO"'"'' cni^'^: 
NGC 1569, NGC 3077). For the remaining galax- 
ies, Haro 3 and II Zw 40, no constraints can be 
made (see §6 for details). "Point-source" ratios 
preferentially bias the ratios downward, ie. to- 
wards cooler temperatures or lower densities, so 
these derived values reflect lower limits to the den- 
sity and temperature if the true source sizes are 
larger than ^ 10 . 

Table 4 summarizes the excitation conditions 
for the eight galaxies. The derived tempera- 
tures and densities of the molecular clouds in 
these dwarf starburst galaxies are different from 
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those found in the Galactic disk and non-starburst 
dwarf galaxies, but similar to those found in other, 
higher metallicity, starburst nuclei (eg. Young and 
Scovillc 1991; Braine and Combes 1992; Wild et 
al. 1992; Devereuxet al. 1994; Aalto et al. 1995; 
Petitpas and Wilson 1998). Sage et al. (1992) 
come to a similar conclusion based on CO (2-1). 

Despite the low metallicity, strong radiation 
fields and weak CO, CO appears not to have 
become optically thin in the low ,J transitions 
of dwarf starbursts. A likely explanation for 
this is that if the CO column density is high 
enough to shield itself against photodissociation 
and hence maintain a detectable amount of CO 
(Nco ~ 10^'^ cm-2; Lequetix et al. 1994), then 
that column density is high enough for the low 
J transitions of CO to be optically thick (Pak 
et al. 1998). Over the region of the starburst 
where CO can survive, CO might be expected to 
be optically thick and resemble the physical condi- 
tions of higher metallicity starbursts, with the dif- 
ference between the two environments (the weak- 
ness of CO in the low metallicity systems) reflect- 
ing, predominately, increased beam dilution due to 
smaller "CO cloud size" (eg. Maloney and Black 
1988; Lequeux et al. 1994; Maloney and Wolfire 
1997). 

The large optical depths, the relatively high 
densities derived from LVG modeling plus the 
small molecular cloud filling factor estimates indi- 
cate that the CO cloud size may be smaller in these 
dwarf starbursts than in spiral starbursts. A simi- 
lar effect is observed for the LMC and SMC (Israel 
et al. 1986; Rubio, Lequeux and Boulanger 1993; 
Mochizuki et al. 1994). However, since none of 
the molecular clouds are resolved here, we cannot 
say whether this is due to photodissociated CO or 
is just a reflection of the gas "traced by CO(3-2)" 
with its larger critical density. 

3.3. Distribution of Warm, Dense Molec- 
ular Gas in Dwarf Starbursts 

Since bright CO(3-2) emission requires warmer 
temperatures and higher densities than C0(1- 
0), it is expected that CO(3-2)/CO(l-0) would 
peak on regions of active star formation. There 
has been some evidence for this in nearby spi- 
rals (Garcia-Burillo, Guelin and Cernicharo 1993; 
Wilson, Walker and Thornley 1997, eg.), as well 
as some of the nearby, less active dwarfs (Petitpas 



and Wilson 1998; Mauersberger et al. 1999). 

For the two galaxies in the sample with spatial 
information, He 2-10 and NGC 5253, we investi- 
gate the spatial variation in the CO(3-2)/CO(l- 
0) ratio. In these two galaxies, CO(3-2)/CO(l-0) 
peaks toward the starburst (Table 5). In the case 
of the fully mapped galaxy. He 2-10, the C0(3- 
2) emission is more strongly peaked towards the 
starburst than is CO(l-O), and the tidal features 
seen in CO (1-0) (peak C) are less pronounced or 
"missing" all together . Baas, Israel and Koorn- 
neef (1994) also find indications that the C0(1- 
0) source size is slightly larger than the C0(3- 
2). Likewise for NGC 5253, the CO(2-1)/CO(1-0) 
value obtained from interferometric data indicates 
that the ratios along the dust lane peak closest to, 
although not at the starburst (Meier and Turner 
1998). The CO (3-2) observations demonstrate a 
similar trend. While the three offset positions 
are not deep enough to detect CO(3-2), the hmits 
show that CO(3-2) is weaker relative to CO(l-O) 
further from the starburst. 

To further test the result that CO (3-2) is 
stronger in the central starburst regions of the 
sample, we compared the observed line ratios 
with the physical size covered on the galaxy. The 
findings are plotted in Figure 5. Note that the 
beamsizes of He 2-10 and NGC 5253 for the mea- 
sured line ratios are larger than 22 and reflect the 
CO(l-O) beam size of 40 and 44 , respectively. 
There is a weak trend (correlation coefficient of r = 
-0.83) that the line ratio decreases as the distance 
of the starburst galaxy increases, which implies 
that including gas in the beam from farther out 
in the galaxy lowers the ratio. (Haro 3 appears 
to be somewhat of an exception, but the line ra- 
tio of this galaxy appears "contaminated" by an 
optically thin velocity component [§6].) This is 
consistent with the notion that the highest tem- 
peratures and densities are concentrated towards 
the starburst, with lower temperature and density 
gas over the rest of the galaxy. We conclude that 
this is a general feature among dwarf starbursts, 
and the properties derived using the higher J tran- 
sitions of CO probably better represent the region 
directly associated with the star formation. 
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4. Masses, Column Densities cind Conver- 
sion Factors 

Estimating molecular gas masses in dwarf 
galaxies can be a particularly tricky business. 
Dwarf galaxies seem to have a higher standard 
conversion factor {Xco, where X = N{H2)/ Ico) 
than do large metal-rich spirals (Israel 1986; Mal- 
oney and Black 1988; Eckart et al. 1990; Verter 
and Hodge 1995; Wilson 1995; Arimoto, So- 
fue and Tsujimoto 1996). Several methods are 
used to obtain molecular mass estimates. One is 
to assume that the molecular clouds are in virial 
equilibrium, and that the linewidth is representa- 
tive of the cloud mass. However, with the large 
beamsizes of single-dish observations, galaxy ro- 
tation contributes significantly to the linewidth. 
In fact, it is more likely that the mass measured 
in this fashion represents an estimate of the dy- 
namical mass over the beam. A second method 
for estimating molecular mass is to assume CO is 
optically thin and "count molecules", using (in- 
cluding He; eg. Rohlfs and Wilson 1996, page 
191): 

M..„ = 173(M,, (^1^) 

[e^ - 1 ) K km s-i 

where the value is scaled relative to the Galac- 
tic [CO/H2] abundance ratio of 8.5 x 10^"'' (Fr- 
erking, Langer and Wilson 1982). This method 
underestimates the total molecular mass, both be- 
cause the CO emission is optically thick (§3.1), 
and the CO abundances, [CO/H2], are probably 
lower than the Galactic value for most of the dwarf 
galaxies. Thirdly, molecular masses can be esti- 
mated using a metallicity scaled conversion factor. 
We use (Wilson et al. 1988; Petitpas and Wilson 
1998): 



M^oi = 1-23 X 10^ (Mq) 




/ 115 GHz V 2 Sco 

V I' ) R ' 

where Xco^ai = 2-3 x lO^" cm-'^{Kkms-^y^ 
(Strong et al. 1988), Sco is the CO (3-2) flux 
in Jy km s~^ (where a conversion factor of 47.3 



Jy/K is assumed), and R is the CO(3-2)/CO(l-0) 
line ratio. The scaling for the conversion factor is 
quite uncertain due to the few data points, but we 
estimate it based on the relationship derived from 
single-dish data (Arimoto, Sofue and Tsujimoto 
1996). 

Collected in Table 6 are the masses derived for 
each galaxy by the three different methods. We 
also compare the derived masses with the dust 
mass estimated from the 100 /im IRAS fluxes 
(Thronson and Telesco 1986; Melisse and Israel 
1994). The virial theorem gives a much larger 
mass than any of the other methods, consistent 
with it tracing dynamical mass. It should give a 
firm upper limit to the molecular mass present. 
The mass estimates obtained by assuming the 
CO (3-2) emission is optically thin are more than 
100 times lower than what is estimated from the 
conversion factor method. The Mthin estimate 
gives a firm lower limit to the amount of molecular 
gas. The conversion factor method is intermedi- 
ate between the other two. Therefore, even with 
its large uncertainty, the best estimate is probably 
still the conversion factor based estimate. Since 
none of the molecular clouds are resolved spatially 
or in velocity, a detailed investigation of the Xco 
must await higher resolution. The CO(3-2) data 
does, however, appear consistent with the Ico 
vs. metallicity relationship found from CO (1-0) 
for dwarf and spiral galaxies (Figure 3 Verter and 
Hodge 1995; Wilson 1995; Arimoto, Sofue and 
Tsujimoto 1996). 

5. CO and Dust: Preferential CO Deple- 
tion? 

A correlation between CO and IR luminosity 
for large starbursts has been known for some time 
(see the Young and Scoville 1991, review and there 
references therein). In Figure 4, the infrared lumi- 
nosity of this sample plus the Mauersberger et al. 
(1999) sample of non-dwarf nuclei are plotted ver- 
sus their CO (3-2) luminosity. The CO (3-2) data 
for the dwarf starbursts also show a tight correla- 
tion between Lco32 and L/ij, with a relation of: 

logiLiR) = (6.31 ± 0.5) + (0.46 ± Om)log{Lco32) 

(correlation coefficient of r=-0.94). This can be 
compared to the correlation we calculate using the 
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Mauersberger et al. (1999) dataset: 



log{LiR) = 3.20 + 0.93log{Lco32) 

In dwarf starbursts, CO emission appears weaker 
relative to the infrared emission than in tlic liigh 
metallicity galaxies. The slope of this relationship 
is substantially steeper than what is found from 
nearby, high metallicity, starburst spirals. The 
most metal-rich galaxies in this sample follow the 
relationship found for non-dwarfs, but the lower 
metallicity sources are weaker in Lco32 relative to 
dust than would be extrapolated from non-dwarfs. 
While the sample of dwarf galaxies is too small to 
prove this fact conclusively, it indicates that there 
may be a break in the slope of Lco32 vs. L//j at 
the faint, low metallicity end. Taylor, Kobulnicky 
and Skillman (1998) also suggest a rapid fall off 
in Leo at metallicities of about 1/10 Zq based on 
a different line of reasoning. Further observations 
are key to investigating the validity of this trend. 

The weakening of T^co relative to L/j? is not 
necessarily unexpected because of the strong radi- 
ation field and low metallicity of these dwarf star- 
bursts. These simultaneously lead to higher dust 
temperatures and hence L//f and to increased CO 
photodissociation. Dust survives longer than CO 
in strong radiation fields, so CO should disappear 
faster than dust. 

6. Notes on the Individual Galaxies 

He 2-10: In Figure 6, the channel maps of He 
2-10 are presented. A 14-point map of He 2-10 was 
made on a 10 grid in RA and Dec. around the 
brighter optical starburst component (Vacca and 
Conti 1992). The (Aa, A(5) = 0,0 position of He 
2-10 was reobserved on February 25 to confirm its 
line strength. Pointings were extended towards 
the southeastern portion of the galaxy where a 
bar-like tidal feature is seen in the CO (1-0) in- 
terferometer map (Kobulnicky ct al. 1995). We 
detect emission in every position except (20 ,0). 
The peak Tmb occurs at the central position and 
has a value of 0.27 K. This value is about 80% of 
the value Baas, Israel and Koornneef (1994) find 
for CO (3-2) after we convolve to the same beam- 
size. The CO(3-2) emission at this position has a 
^FWHM ~ 57 km s~^ centered at 849 km s~^, 
with detectable emission spanning a velocity range 
of 810 - 890 km s'^ (LSR). 




5 ■ 

8 8.5 9 9.5 10 10.5 11 11.5 

log(LJ 

Fig. 5.— The luminosity of CO (3-2) in units 

of K km s^^pc^^ plotted versus the infrared lu- 
minosity in units of L©, for the dwarf galaxies 
(bold pluses) and for a sample of non-dwarf galac- 
tic nuclei (Mauersberger et al. 1998; crosses). 
The dashed line is the least-squares fit to the 
dwarf galaxy sample and the solid line is the least- 
squares fit to the spiral nuclei sample. 
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Fig. 6. — The channel maps of He 2-10. Panels are 
displayed in intervals of 8 km s"'^. The Vlsr of 
each channel is displayed in the upper right hand 
corner of each frame. Contours are 1, 2 8 times 
the 2f7 temperature of 0.035 K. The asterisk marks 
the location of the optical starburst and center 
position (0,0). Each additional pointing is marked 
by a cross. The feature denoted by Kobulnicky 
et al. (1995) as "C", is also marked at its correct 
location and velocity, and the location and the two 
velocities contributing to the "detached" feature is 
marked with a "D" (Kobulnicky et al. 1995). 
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Fig. 7. — The integrated intensity map of He 2-10. 
The integrated intensity map was mapped from 
only the channels displayed in Figure 2, (773 to 
924 km s~^). Emission below l.Scr in each chan- 
nel map were cut. Contours are 1, 2 10 times 
the 2(7 value of 1.6 K km s~^. Again, the aster- 
isk marks the location of the optical starburst and 
center pointing (0,0) and each additional pointing 
is marked by a cross. The map is uncertain at 
the (20 ,-20 ) position because there is no obser- 
vation there. To represent this fact the contours 
are dashed. 



NSC Cd/Red 




Fig. 8.— The spectra of NGC 5253 overlaid 
on the OVRO millimeter interferometer C0(1- 
0)/Red image of Turner, Beck & Hurt (1997) 
(their Figure 2). Each observation is displayed 
with a line running from the spectrum to the loca- 
tion at which it was obtained. The vertical axis is 
main-beam temperature and ranges from -0.065 to 
0.065 K. The horizontal axis is velocity and ranges 
from 180 to 620 km s~^. The spectra are hanning 
smoothed to a velocity resolution of 8 km s~^. For 
the central position the new CO (2-1) spectrum is 
included imderneath the CO (3-2) spectra. For the 
detections the best fit gaussian is displayed. 



10 



Wc do not detect any obvious emission to a level 
of 40 mK associated with the kinematically dis- 
tinct feature labeled by Kobulnicky et al. (1995) 
with the letter "C" . There is some hint of a south- 
eastern extension seen at 853 km s~^ and 861 km 
s~^, but this is separated by better than two chan- 
nels (16 km s~^) from their velocity and is not 
likely to be the same source. On the other hand, 
we do detect emission that matches both in posi- 
tion and velocity space with the weak detached 
feature seen in the Kobulnicky et al. (1995) 
interferometric CO(l-O) map. Kobulnicky et al. 
(1995) interpret this detached feature question- 
ably as noise, but since we see indications of emis- 
sion in CO (3-2) at the same location and veloc- 
ities, this feature may be real. This feature is 
weak, with a Tmb ~ 0.08 K, and centered at 
a = 8''34™08.*65 and 5 = - 26°13'56."5. 

Figure 7 shows the integrated CO (3-2) inten- 
sity map of the nuclear region of He 2-10. The 
emission from the central peak is more dominant 
than in CO(l-O). Within pointing uncertainties, 
the peak of the integrated intensity coincides with 
the position of the brighter optical starburst com- 
ponent. Of the three extensions seen, one extends 
to the northeast along a p. a. ~ 50°, while the 
other two extend to the southeast. We interpret 
the two extensions to the SE as one feature along a 
p. a. of 130° and attribute the apparent split into 
two features to the lack of a pointing at (20 ,- 
20 ). This extension matches what is observed 
in the interferometric CO(l-O) and C0(2-l) maps 
and is near the site of an offset star cluster (Kob- 
ulnicky et al. 1995; Beck and Kovo 1998; Meier 
and Turner 1998). 

The centrally peaked component has a FWHM 
of 26 . When deconvolved from the beam, we es- 
timate that the source size is 13 ± 3 (560±110 
pc at 9 Mpc). Since the observations of Baas, 
Israel and Koornneef (1994) are at a different res- 
olution than ours, we can compare the observed 
main-beam temperatures to get a second estimate 
of the source size. Again, 13 is obtained. Because 
both methods agree, and are consistent with the 
size estimated from the high resolution interferom- 
eter C0(2-l) map (Meier and Turner 1998), we 
conclude that this is a good representation of the 
CO(3-2) source size. 

Assuming a source size of 13 for He 2-10, we 
obtain CO(3-2)/CO(2-l), CO(3-2)/CO(l-0), and 



CO(2-1)/CO(1-0) line ratios of 1.0±0.2, 0.61±0.1 
and 0.59±0.1 respectively. These are slightly 
lower than the previous results of Baas, Israel and 
Koornneef (1994). CO(3-2)/CO(2-l) > 1.0 com- 
bined with a high ^^CO/^^CO isotopic ratio led 
them to conclude that the CO emission is opti- 
cally thin. With the advantage of a fully sam- 
pled map and CO(3-2) observations with beam 
matched to their J=2-l observations, these obser- 
vations do not necessarily require optically thin 
gas. LVG modeling yields a best fit solution of 
> 10^-^ cm~^ and Tk ~ 5 - 10 K, however the 

on the solution is high. These gas temperatures 
are slightly cooler than the ~20 K estimated from 
LTE modeling. The 12CO(1-0)/13CO(1-0) line ra- 
tio was also modeled. A model with [^^CO/^^CO] 
= 40 and dv/dr = 1 km s~^pc~^ and a ^'^CO de- 
pleted model with \^^CO/'^^CG\ = 120 were run 
(Figure 2). Both CO abundance models arc in- 
consistent with the solution found from the ^^CO 
transitions. The ^^CO line in He 2-10 is too weak 
(i2CO(l-0)/i3CO(l-0) < 22) relative to what is 
predicted from the ^^CO lines, as Baas, Israel 
and Koornneef (1994) noted. In order to have 
i2CO(l-0)/"CO(l-0) - 22, the i^co lines must 
be optically thin, unless the abundance is severely 
depleted. However nowhere over the CO (3-2) so- 
lution space is the modeled -'^^CO optically thin. 
Therefore, given the high on the ^^CO solution 
and that the ^^CO observations are not consis- 
tent with the simple LVG modeling, the solution 
should be considered uncertain. ^^CO is proba- 
bly preferentially photodissociated or that source 
sizes arc much smaller in -'^^CO than -'^^CO. As far 
as the high isotopic ratios obtained are concerned 
(Baas, Israel and Koornneef 1994; Kobulnicky et 
al. 1995), the optically thinner '^"^CO lines are 
less capable of self- shielding than ^^CO, and are 
likely to be preferentially photo-dissociated by the 
intense radiation fields. So the abimdancc and 
spatial extent of ^"^CO may be greatly diminished 
with respect to ^^CO (eg. van Dishoeck and Black 
1988; Lequeux et al. 1994; Warin, Benayoun and 
Viala 1996). 

NGC 5253: The four pointings for NGC 5253 
were chosen based on the CO(l-O) interferometric 
map made at Owens Valley Millimeter Observa- 
tory (OVRO) (Turner, Beck and Hurt 1997). We 
observed two pointings associated with the cen- 
tral dust lane, and two toward a possible molecu- 
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lar cloud at the edge of the OVRO primary beam. 
In Figure 8, the four CO(3-2) and one C0(2-l) 
spectra are overlaid on the OVRO CO(l-O) map 
(Figure 2 of Turner, Beck and Hurt 1997). The 
measured CO (3-2) intensity is about 5 - 6 a, and 
both the velocity centroid (Vo ~ 420 km s~^) and 
the line width (AVfwhm ^ 45 km s"-'^) match 
published single-dish CO(l-O) detections (Wiklind 
and Henkel 1989; Taylor, Kobulnicky and Skill- 
man 1998). For the three other spectra there are 
no detections to 3a. LVG modeling yields den- 
sities of > 10^-^ cm~^, and T^ > 40 K, 
with the best fit solutions obtained for the lower 
densities. 

NGC 1569: We detect CO(3-2) in this galaxy 
with Ico=2-2 K km s~^ {Tmb = 63 mK) towards 
the HII region C (Waller 1991). This corresponds 
to the central position observed by Greve et al. 
(1996) in CO(l-O) and C0(2-l). The line width 
of CO(3-2) is shghtly larger (32.3 km s'^) than 
that of the CO(l-O) observations (23.6 km s'^). 
The line center is shifted to slightly higher veloci- 
ties than the CO(l-O) (-69 relative to -81 km s'^ 
for CO(l-O)). The C0(2-l) hne center is inter- 
mediate between the two, 78.4 km s~^. There- 
fore, while the velocity centroids are roughly con- 
sistent, there is a marginal trend that the emission 
from the higher J transitions tend towards increas- 
ing velocity. Comparing these higher velocities 
with the velocity field obtained from HI (Reakes 
1980; Stil and Israel 1998) indicates that C0(3- 
2) peaks slightly closer to the super-star cluster, 
A, than does CO(l-O), and is probably associated 
with CMC 2 of Taylor et al. (1999). 

Interestingly, there is a secondary peak seen at 
-140 km s^^. This peak is at the 4cr level, and the 
line profile is similar to the main peak. HI is also 
present at this velocity, giving more evidence this 
feature may be real. We have inspected the C0(1- 
0) and C0(2-l) spectra of Yoimg, Gallagher and 
Hunter (1984), Taylor, Kobulnicky and Skillman 
(1998) and Greve et al. (1996) for similar fea- 
tures. Unfortunately, the spectra of Greve et al. 
(1996) does not extend that far in velocity, and 
while the Young, Gallagher and Hunter (1984) 
spectrum shows indications of a feature at this ve- 
locity, it does not seem to be confirmed by the 
more sensitive observations of Taylor, Kobulnicky 
and Skillman (1998). Therefore, until further ob- 
servations can be obtained, we assume it is not 



real. Since both the CO(3-2)/CO(l-0) and the 
CO(2-1)/CO(1-0) line ratios are ~ 1, a large re- 
gion of parameter space is acceptable to the LVG 
models. The best fits represent an arch running 
from Tk ~ 100 K and — 10^'^ cm~^ down 
to Tk ~ 40 K and ~ 10^-^ cm~^ and back 
over to Tk 100 K and njy^ ~ 10^-^ cm,-^. 

NGC 3011: Observations were centered on 
the location of the interferometer CO (2-1) peak 
(Meier and Turner 1998). This is equivalent to 
the (-10,0) position of Becker, Schilke and Henkel 
(1989). The CO(3-2) line width of 51 km is 
slightly wider than observed by Becker, Schilke 
and Henkel (1989) in CO(l-O) (34 km s"^), while 
the centroid is consistent with what is observed 
in the lower transitions. Based on the C0(2- 
1)/CO(1-0) hne ratio discussed in Becker, Schilke 
and Henkel (1989), we estimate a CO(3-2)/CO(l- 
0) line ratio of 1.1. Becker, Schilke and Henkel 
(1989) obtain a CO(2-1)/CO(1-0) line ratio of 
0.82. For NGC 3077, the best fit LVG solution 
corresponds to T^ ~ 30 K and n//^ > 10^'^ cm~^, 
with an acceptable range of solutions covering 
Tk > 20 K and nn^ > 10^ cm,-^. 

Haro 2: We observed CO(3-2) at the same 
position as CO(l-O) and C0(2-l) of Sage et al. 
(1992). The fitted line width is 69 km s'^ but 
the line appears distinctly non-gaussian. The SNR 
of the spectrum is to low to warrant any multi- 
component fit, but the spectrum is consistent with 
two roughly equal intensity gaussians separated by 
about 40 km s~^. For Haro 2, the best fit solution 
obtained from LVG modeling is n//^ > lO"' cm~^ 
and Tfc ~ 5 K, but the are high, similar to 
what is fomid for He2-10. However, the low ra- 
tios obtained for both transition require that T^ 
be low and riH^ be high. 

Haro 3 (NGC 3353): C0(3-2) observations are 
centered on the starburst region B (Steel et al. 
1996). The fitted gaussian yields a line width of 
53 km s""'^ which is slightly narrower than C0(1- 

0) (A^;io = 64 km s"^; Sage et al. 1992). The 
CO (3-2) spectrum shows a narrow brighter com- 
ponent at 960 km s"'^ {vlsr = 950 km s"-'^). The 
same peak is seen in both the CO (1-0) and C0(2- 

1) spectra of Sage et al. (1992), and matches the 
velocity of Ha at the starburst region, Bl, of Steel 
et al. (1996). Therefore, we conclude this feature 
is real. Inspecting the spectra of Sage et al. (1992) 
a main-beam temperature ratio of 1.3 is estimated 
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for this component. The large CO(3-2)/CO(l-0) 
line ratio is inconsistent with any region of param- 
eter space indicated by the CO(2-1)/CO(1-0) hne 
ratio, so no LVG solutions is foimd. This compo- 
nent is potentially optically thin and may be caus- 
ing an overestimate of the global CO(3-2)/CO(l- 
0) ratio in Haro 3. 

// Zw 40 : The beam was centered on the star- 
formation peak and is within an arcsecond of the 
observed C0(2-l) and CO(l-O) position of Sage 
et al. (1992). The line is not detected. Our 
limit for the CO (3-2) main-beam temperature is 
7 mK. There appears to be no hint of a line at 
the 1(7 level in either of the previously published 
CO(l-O) velocity centroids, ~ 770 km s~^ (Sage et 
al. 1992) or ~ 850 km (Tacconi and Young 
1987). LVG modeling does not provide a useful 
constraint on the gas properties because CO(3-2) 
is not detected. The low CO(3-2)/CO(l-0) ratio 
implies that the molecular gas must not be very 
warm and dense. 

Mrk 86 (NGC 2537): We do not detect any 
CO(3-2) emission to a limit of 11 mK. Sage et 
al. (1992) detected weak CO(l-O) and C0(2-l) 
at a Vlsr of 460 km s~^. There is a hint of a very 
weak line at Ylsr = 428 km s"'^, but only at 1.5cr. 
The low CO(3-2)/CO(l-0) ratio limit and the low 
CO(2-1)/CO(1-0) ratio imply LVG fits quite sim- 
ilar to those of Haro 2. 

7. Conclusions 

We have observed CO (3-2) in a sample of eight 
dwarf galaxies with very intense recent star for- 
mation. The galaxies range in metallicity from 1 
Z© to ~0.1 Zq. C0(3-2) was detected in six of 
the eight galaxies. The galaxies are all quite weak 
(Tmb ~ 50 mK), with the exception of He 2-10, 
which was mapped. The two galaxies not detected 
(II Zw 40 and Mrk 86) are two of the lowest metal- 
licity galaxies. 

The CO(3-2)/CO(l-0) line ratio ranges from 
0.37 - 1.1 for the six detected galaxies. With the 
possible exception of one velocity component in 
Haro 3, all the galaxies appear dominated by op- 
tically thick CO emission, even though they have 
low metallicity and strong radiation fields. This is 
consistent with the fact that the only place that 
CO can survive intense radiation fields is where it 
is self-shielded, and thus optically thick. 



The error weighted mean CO(3-2) /CO(l-O) Hne 
ratio is 0.60±0.06. This is virtually identical to 
what is seen in high metallicity star bursts. While 
CO(3-2) is very weak on an absolute scale, it is 
as bright relative to CO(l-O) as in high metallic- 
ity starbursts. In general, in terms of the physical 
properties of the molecular gas, dwarf starbursts 
appear not to be particularly unusual examples of 
the starburst family. They appear to be similar to 
their high metallicity counterparts, with the ex- 
ception that CO has a smaller filling factor and 
hence weaker intensities. 

LVG models from the CO line ratios imply that 
for He 2-10, Haro 2 and Mrk 86, the molecular 
gas traced by CO is cool (T^ ~10 K) and dense 
(nH,< 10^ cm-3). For NGC 5253, the molecular 
gas appears warm (T^ > 40 K) and moderately 
dense (n//^ ~ 10^'^"^ cm~^). The molecular gas 
in NGC 1569 and NGC 3077 are also warm but 
somewhat more dense {nH2 > 10^ cm~^). 

For the two galaxies that have been observed 
with multiple pointings. He 2-10 and NGC 5253, 
we find that the CO(3-2)/CO(l-0) line ratio de- 
creases away from the starburst. The entire sam- 
ple shows a trend where the ratio is lower for 
the galaxies where the beam includes more "non- 
starburst" gas (ie. covers larger fractions of the 
galaxy). This is interpreted as due to a drop in 
gas temperature and density away from the star- 
burst. 

The combination of low filling factors derived 

from the CO (3-2) emission, the fact that CO (3-2) 
systematically remains optically thick in spite of 
low metallicities and strong radiation fields, and 
a trend found between Leo and Ljr implies the 
lower metallicity galaxies are more strongly de- 
pleted in CO relative to dust than higher metal- 
licity galaxies. 

We gratefully thank Tom Phillips, Dominic 
Benford and Atilla Kovacs for useful discussions 
and help in making the observations. We thank 
Gene Serabyn for obtaining the CO (2-1) spectrum 
of NGC 5253. We also thank Chip Kobulnicky 
for useful discussions and providing us with the 
OVRO CO(l-O) data for He 2-10. We also thank 
the referee, Dr. Christine Wilson for many helpful 
comments. This work was supported in part by 
NSF grant AST9417968 and by the Sackler Cen- 
ter for Astronomy at Tel Aviv University. This 



13 



research has made use of the NASA/IPAC extra- 
galactic database (NED). 

REFERENCES 

Aalto, S., Booth, R. S., Black, J. H. & Johansson, 
L. E. B. 1995, A&A, 300, 369 

Arimoto, N., Sofue, Y. & Tsujimoto, T. 1996, 
PASJ, 48, 275 

Arnault, P., Casoh, F., Combes, F., & Kunth, D. 
1988, A&A, 205, 41 

Arp, H. & Sandage, A. 1985, AJ, 90, 1163 

Baas, F., Israel, F. P. & Koornneef, J. 1994, A&A, 
284, 403 

Beck, S. C. & Kovo, O. 1998, AJ, 117, 190 

Becker, R., Schilke, P. & Henkel, C. 1989, A&A, 
L211, 19 

Braine, J. & Combes, F. 1992, A&A, 264, 433 

Conti, P. S. 1991, ApJ, 377, 115 

Conti, P. S. & Vacca, W. D. 1994, ApJ, 473, L97 

Cottrell, G. A. 1976, MNRAS, 174, 455 

De Jong, T., Chu, S.-I. & Dalgarno, A. 1975, ApJ, 
199, 69 

Devereux, N., Yoshiaki, T., Sanders, D. B., Nakai, 
N. & Young, J. S. 1994, AJ, 107, 2006 

Eckart, A., Downes, D., Gcnzcl, R., Harris, A. I., 
Jaffe, D. T.,& Wild, W. 1990, ApJ, 348, 434 

Elmegreen, B. G. 1989, ApJ, 338, 178 

Frerking, M. A., Langer, W. D. & Wilson, R. W. 

1982, ApJ, 262, 59 

Gallagher, J. S. & Hunter, D. A. 1984, ARA&A, 
22, 37 

Garcia-Burillo, S., Guelin, M. & Cernicharo, J. 
1993, A&A, 274, 123 

Garnett, D. R. 1990, ApJ, 363, 142 

Goldreich, P. & Kwan, J. 1974, ApJ, 189, 441 

Gonzalez Delgado, R. M., Leitherer, C, Heckman, 
T. & Cervino, M. 1997, ApJ, 483, 705 



Grave, A., Becker, R., Johansson, L. E. B. & McK- 
eith, C. D. 1996, A&A, 312, 391 

Heckman, T. M. 1980, A&A, 87, 142 

Hodge, P. W. 1971, ARA&A, 9, 35 

Hurt, R. L., Turner, J. L., Ho, P. T. P. & Martin, 

R. N. 1993, ApJ, 404, 602 

Israel, F. P. 1986, A&A, 168, 369 

Israel, F. P. 1988, A&A, 194, 24 

Israel, F. P., de Bruyn, A. G. 1988, A&A, 198, 109 

Israel, F. P., De Gwaauw, Th., Van de Stadt, H. 
& De Vries, C. P. 1986, ApJ, 303, 186 

Israel, F. P., & van Driel, W. 1990, A&A, 236, 323 

Klein, U., Weiland, H., & Brinks, E. 1991, A&A, 
246, 323 

Kobulnicky, H. A., Dickey, J. M., Sargent, A. I., 
Hogg, D. E. & Conti, P. S. 1995, AJ, 110, 116 

Kobulnicky, & H. A., Johnson, K. 1999, ApJ, 527, 
154 

Kobulnicky, H. A., Kennicutt, R. C, & Pizagno, 
J. L. 1998, ApJ, 514, 544 

Kobulnicky, H. A., & SkiUman, E. D. 1995, ApJ, 

L454, 121 

Kobulnicky, H. A., & SkiUman, E. D. 1997, ApJ, 

489, 636 

Kunth, D., & Jobert, M. 1985, A&A, 142, 411 

Lequeux, J., Le Bourlot, J., Pineau Des Forets, 
G., Roueff, E., Boulanger, F. & Rubio, M. 1994, 
A&A, 292, 371 

MacLaren, I., Richardson, K. M. & Wolfendale, 
A. W. 1988, ApJ, 333, 821 

Maloney, P. & Black, J. H. 1988, ApJ, 325, 389 

Maloney, P., & Wolfirc, M. G. 1997, in lAU Symp 
170: CO: Twenty-Five Years of Millimeter- 
Wave Spectroscopy, ed. W. B. Latter et al. 
(Dordrecht:Kluwer), 299 

Marconi, G., Matteucci, F., & Tosi, M. 1994, MN- 
RAS, 270, 35 



14 



Maucrsbcrgcr, R., Henkel, C, Walsh, W. & 
Schulz, A. 1999, A&A, 341, 256 

Meier, D. S., & Turner, J. L. 1998,in: K.A. van 
der Hucht, G. Koenigsbcrger & P.R.J. Eenens 
(eds.), Wolf-Rayct Phenomena in Massive Stars 
and Starburst Galaxies, Proc. lAU Symp. No. 
193 (San Francisco: ASP), 746 

Meier, D. S. & Turner, J. L. 2000, in prep. 

Meier, D. S., Turner, J. L., & Hurt, R. L. 2000, 
ApJ, 531, 200. 

Melisse, J. P. M. & Israel, F. P. 1994, A&AS, 103, 
391 

Meurer, G. R., Heckman, T. M., Leitherer, C., 
Kinney, A., Robert, C. & Garnett, D. R. 1995, 

AJ, 110, 2665 

Mochizuki, K., et al. 1994, ApJ, 430, L37 

Niklas, S., Klein, U., Braine, J., & Wielebinski, R. 
1995, A&AS, 114, 21 

Pak, S., Jaffee, D. T., van Dishoeck, E. P., Jo- 
hansson, L. E. B., & Booth, R. S. 1998, ApJ, 
498, 735 

Pctitpas, G. R. & Wilson, C. D. 1998, ApJ, 496, 

226 

Reakes, M. 1980, MNRAS, 192, 297 

Rohlfs, K. & Wilson, T. L. 1996, Tools of Radio 
Astronomy, 2nd Ed., (Springer- Verlag:Bcrlin) 

Rubio, M., Lequeux, J. & Boulanger, F. 1993, 
A&A, 271, 9 

Sage, L. J., Salzer, J. J., Loose, H.-H. & Henkel, 
C. 1992, A&A, 265, 19 

Sandage, A. 1994, ApJ, L423, 13 

Sanders, D. B., Tilanus, R. P. J., Scoville, N. Z., 

Wang, Z., & Zhou, S. 1993, in Back to the 
Galaxy, ed. F. Verter (Kluwer, Dordrecht), 21 

Schaerer, D., Contini, T. & Pindao, M 1999, 
A&AS, 136, 35 

Skillman, E. D., Kennicutt, R. C. & Hodge, P. W. 
1989, ApJ, 347, 875 



Steel, S. J., Smith, N., Metcalfe, L., Rabbette, M. 
& McBreen, B. 1996, A&A, 311, 721 

Stil, J. M. & Israel, F. P. 1998, A&A, 337, 64 

Strong et al. 1988, A&A, 207, 1 

Tacconi, L. J. & Young, J. S. 1987, ApJ, 322, 681 

Tammann, G., A., & Sandage, A. 1968, ApJ, 151, 

825 

Taylor, C. L., Hiittcmcistcr, S., Klein, U., & 
Grevc, A. 1999, A&A, 349, 424 

Taylor, C. L., Kobulnicky, H. A. & Skillman, E. 
D. 1998, AJ, 116, 2746 

Thronson, H. A. & Telesco, C. M. 1986, ApJ, 311, 
98 

Thuan, T. X. & Martin, G. E. 1981, ApJ, 247, 823 

Turner, J. L., Beck, S. C. & Hurt R. L. 1997, ApJ, 
L474, 11 

Turner, J. L., Ho, P. T. P., & Beck, S. C. 1998, 
AJ, 116, 1212 

Turner, J. L., Hurt R. L., & Hudson, D. Y. 1993, 
ApJ, L413, 19 

Vacca, W. D. & Conti, P. S. 1992, ApJ, 401, 543 

van Dishoeck, E. F., & Black, J. H. 1988, ApJ, 
334, 771 

Verter, F. & Hodge, P. 1995, ApJ, 446, 616 

Waller, W. H. 1991, ApJ, 370, 144 

Warin, S., Benayoun, J. J., & Viala, Y. P. 1996, 
A&A, 308, 535 

Wiklind, T. & Henkel, C. 1989, A&A, 225, 1 

Wild, W., Harris, A. I., Eckart, A., Genzel, R., 
Graf, U. U., Jackson, J. M., Russell, A. P. G. 
& Stutzki, J. 1992, A&A, 265, 447 

Wilson, C. D. 1995, ApJ, L448, 97 

Wilson, C. D., Howe, J. E, & Balogh, M. L. 1999, 
ApJ, 517, 174 

Wilson, C. D., Scoville, N. Z., Freedman, W. L., 
Madore, B. F., & Sanders, D. B. 1988, ApJ, 
333, 611 



15 



Wilson, C. D., Walker, C. E. & Thornley, M. D. 
1997, ApJ, 483, 210 

Wolfire, M. G., Tielens, A. G. G. M. & HoUenbach, 
D. J. 1990, ApJ, 358, 116 

Young, J. S., Gallagher, J. S. & Hunter, D. A. 
1984, ApJ, 276, 476 

Young, J. S. & Scoville, N. Z. 1991, ARA&A, 29, 
581 



This 2-column preprint was prepared with the AAS lAlJ^X 
macros v5.0. 



16 



Table 1 



The Galaxy Sample 


Galaxy 


RA 


Dist. 


IRAS eO/lOOyum 




logfMwr) 


References 




DEC 


(Mpc) 


(Jy) 


V W / 


V w / 










J- D 


i.\±jj 








(B1950) 




(K) 


(Mo) 






He 2-10 


08:34:07.3 


9.0 


24.0/26.4 


9.43 


8.53 


1,9,10,11 




-26:14:05.9 


17.4 


43 


5.47 


8.93 




NGC 5253 


13:37:05.8 


4.1 


30.9/29.0 


8.98 


8.3 


12,15,19,22,23 




-31:23:19.0 


17.2 


46 


4.73 


8.10 




NGC 1569 


04:26:02.0 


2.2 


46.3/50.7 


8.65 


7.93 


5,6,7,13,22 




64:44:31.0 


16.9 


43 


4.52 


8.19 




NGC 3077 


09:59:20.0 


3.25 


14.7/26.9 


8.70 


9.00 


2,4,16,17,21 




68:58:30.0 


16.1 


34 


4.99 


9.02 




Haro 2 


10:29:22.7 


20.3 


4.7/5.3 


9.60 


8.68 


1,14,18,22 




54:39:24.0 


18.4 


42 


5.51 


8.4 




Haro 3 


10:42:16.5 


13.7 


5.1/6.4 


9.30 


8.79 


8,18,20,22 




56:13:23.0 


17.5 


41 


5.16 


8.4 




II Zw 40 


05:53:04.9 


9.2 


6.5/5.7 


9.04 


8.3 


3,8,18,22 




03:23:06.0 


16.2 


48 


4.66 


8.15 




Mrk 86 


08:09:42.8 


6.3 


3.2/6.3 


8.45 


8.28 


8,18,22 




46:08:33.0 


16.6 


33 


5.00 


8.05^^ 





"^using R(60/100)=0.6-^(ei3*5/^^-l)/ (e'^'^^/'^^-l), (Thronson & Telesco 1986) 
''using MD=5D'^Fjy{e^^'^/'^^-l), (Th ronson & Telesco 1986) 

'^Extrapolated from the absolute magnitude based on the relationship of Skillman, Kennicutt & 
Hodge 1989 

References.— (l)Baas et al. 1994; (2)Cottrell 1976; (3)Garnettl990; (4)Heckman 1980; (5)Is- 
rael 1988; (6)Israel & de Bruyn 1988; (7)Israel & van Driel 1990; (8)Klcin et al. 1991; (9)Kobulnicky 
et al. 1995; (lO)Kobulmcky & Johnston 1999; (ll)Kobulnicky et al. 1998; (12)Kobulnicky & Skill- 
man 1995; (13)Kobulmcky & Skillman 1997; (14) Kunth & Jobert 1985; (15)Marconi et al. 1994; 
(16)Melisse & Israel 1994; (17)Niklas et al. 1995; (18)Sage et al. 1992; (19)Sandage 1994; (20)Steele 
et al. 1996; (21)Tammann & Sandage 1968; (22)Thronson & Telesco 1986; (23)Turner et al. 1998 
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Table 2 
Results 



Galaxy 


Transition 


Offset 
(Aa", as") 


(km s~^) 


AV 
(km s-i) 


K 


K km 


He 2-10 


3-2 


(0,0) 


849±1.0 


56.6±2.3 


0.27±0.018 


16.6±0.6 






(0,-10) 


846±1.6 


63.8±3.7 


0.20±0.018 


13.3±0.7 






(-10,0) 


854±1.5 


58.8±3.6 


0.19±0.02 


11.9±0.6 






(0,10) 


850±2.0 


63.4±4.8 


0.15±0.022 


10.3±0.7 






(10,0) 


846±1.4 


50.0±3.4 


0.24±0.022 


12.6±0.7 






(10,10) 


849±3.0 


61.7±7.0 


0.12±0.016 


7.77±0.8 






(10,-10) 


833±1.8 


37.2±4.4 


0.19±0.023 


7.47±0.7 






(10,-20) 


872±14.4 


140±25 


0.047±0.024 


7.00±1.4 






(20,10) 


826±3.3 


22.5±7.1 


0.090±0.022 


2.16±0.6 






(0,-20) 


845±5.6 


33.0±4.4 


0.072±0.030 


2.53±0.8 






(-10,-10) 


851±2.3 


59.7±6.5 


0.142±0.020 


9.06±0.8 






(20,0) 






<0.053 








(-10,10) 


855±3.0 


73.7±7.7 


0.094±0.019 


7.4±0.6 






(20,-10) 


825±3.9 


26.7±10.8 


0.103±0.022 


2.91±0.7 


NGC 5253 


3-2 


(0,0) 


421±4.2 


45.4±8.7 


0.070±0.021 


3.39±0.6 






no -10"! 






<0 07"^ 








(20,20) 






<0.069 








(15,25) 






<0.075 






2-1 


(0,0) 


407±2.9 


71.6±9.0 


0.034±0.006 


2.58±0.3 


NGC 1569 


3-2 


(0,0) 


-68.4±2.0 


32.3±4.0 


0.063±0.01 


2.2±0.3 


NGC 3077 


3-2 


(0,0) 


7.80±2.9 


50.9±7.0 


0.10±0.02 


5.5±0.6 


Haro 2 


3-2 


(0,0) 


1422±4.6 


68.8±9.8 


0.032±0.01 


2.3±0.3 


Haro 3 


3-2 


(0,0) 


951±4.6 


52.6±14 


0.034±0.01 


1.9±0.4 


II Zw 40 


3-2 


(0,0) 






<0.021 


<0.9 


Mrk 86 


3-2 


(0,0) 






<0.033 


<1.8 



Note. — Offsets are in arcseconds from the central position listed in Table 1. Errors are Icr 
from a gaussian fit, based only on the imcertainty in the spectrum. The noise in the gaussian 
fit for Haro 2 is worse than the actual noise is the spectrum, implying it is not gaussian shaped 
(see text). For non-detections 3cr limits are quoted. The uncertainties in absolute calibration 
were not included. 



18 



Table 3 
CO Line Ratios 



Galaxy 


Method 


3-2/1-0 


3-2/2-1 


2-1/1-0 


He 2-10 


point source (PS) 


0.50±0.1 


0.87±0.2 


0.48±0.07 




iTnifnrm fillinp' Mil 


1.7±0.3 


96±0 2 


1.73±0.3 




deconvolved (D) 


*0.61±0.1 


0.96±0.2 


*0.59±0.1 


NGC 5253 


PS 


0.65±0.2 


0.70±0.2 


0.93±0.2 


(0,0) 


U 


2.6±0.7 


1.3±0.4 


2.0±0.4 




D 


*0.70±0.2 


0.73±0.2 


*0.96±0.25 


NGC 1569 


PS 


*1.00±0.2 


2.5±0.6 


0.42±0.1 




U 


1.10±0.2 


0.88±0.2 


1.10±0.2 




D 






*1.10±0.2 


NGC 3077 


PS 


*0.95±0.2 


2.1±0.7 


0.55±0.1 




U 


1.04±0.25 


0.83±0.2 


1.2±0.3 




D 






*0.82±0.2 


Haro 2 


PS 


*0.34±0.1 


1.1±0.2 


*0.31±0.10 




U 


0.37±0.1 


0.96±0.2 


l.OiO.l 


Haro 3 


PS 


*0.89±0.3 


1.7±0.5 


*0.50±0.09 




U 


0.97±0.3 


0.57±0.2 


1.7±0.3 


II Zw 40 


PS 


<0.39 




0.58±0.16 




u 


<0.43 


<0.24 


1.78±0.32 


Mrk 86 


PS 


<0.41 




0.31±0.11 




u 


<0.45 


<0.52 


1.0±0.36 



Note. — See the text for the description of the method for finding each 
ratio and the estimation of the errors. The slight difference between the 
IRAM 30m (1-0) beam and the CSO (3-2) beam has been accounted for 
in the calculations of the ratios. An asterisk denotes which value of the 
line ratio is used in the LVG modeling. All line ratios are based on a 
22 bcamsizc except He 2-10 and NGC 5253, which arc based on a 40 
beamsize. The data for the 2-1 & 1-0 emission comes from: He 2-10: 
Baas et al. 1994; NGC 5253; Wiklind & Henkel 1989 and this paper; 
NGC 1569: Greve et al. 1996; NGC 3077: Becker et al. 1989; the rest: 
Sage et al. 1992. 
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Table 4 
Molecular Gas Properties 





He 2-10 


NGC5253 


NGC1569 


NGC3077 


Haro2 


Haro3 


IIZw40 


Mrk86 


LTE: 


















Td 


43 


46 


43 


34 


42 


41 


48 


33 


T 

J- ex 


11 


15 


>50 


>50 


6 


>50 


<7 


<7 


fa 


0.02 


0.005 


<0.001 


<0.002 


0.005 


<0.001 


<0.004 


<0.005 




103.4 


10^ 


102.5 


10^-7 


103.6 


103.3 


10^ 


102-4 


LVG: 




















10^ 


30 


>40 


^20 


5 






5 




>103-5b 


~ 10^ 


> 102-5 


>103 


>104 








s/n 
Tk 


>5.6b 


1.0 


~1.5 


>1.5 


20 






20 




>3 X lO^b 


3 X 10* 


>1 X 10^ 


>2 X 10"* 


>5 X 10^ 






>5 X 10^ 



'^Using the "OB association" relationship, < Go >— 1-2 x 10~^ Lm of Wolfire, Tielens & HoUenbach 
1990, with R = A ~ 100 pc which is typical of galactic nuclei. This may be an underestimate of A and 
thus of < Go >, since these dwarfs arc dust poor and are potentially more "porous" to radiation than 
metal rich galaxies (Pctitpas & Wilson 1998). 

''The of this solution are high, also the ^"^CO observations are inconsistent with these values (Figure 
2a; §6), so this solution should be considered uncertain. 
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Table 5 
Off-Center Line Ratios 



Galaxy 


Location 


CO(3-2)/CO(l-0) 


He 2-10 


(20,-10) 


0.35 ±0.1 


(D) 


(20,10) 


0.30 ±0.1 


(C) 


(10,-10) 


0.35 ±0.1 


NGC 5253 


(10,-10) 


< 0.7 




(20,20) 


< 0.6 




(15,25) 


< 1.2 



Note. — Ratios are based on the CO(l-O) 
interferometer maps at 22 resolution, and 
have been corrected for resolved out flux (He 
2-10: 82%; Kobulnicky et al. 1995; NGC 5253: 
50%; Turner et al. 1997). Uncertainties for the 
CO (3-2) data are given in Table 3, while the 
interferometric data is assumed to have an un- 
certainty of 20%. In the case of NGC 5253 the 
ratio is of peak Tmb since these three C0(3- 
2) pointings are non-detections: quoted values 
are 3a upper limits. 



Table 6 
Gas Masses 



Galaxy 


T 

ex 

(K) 


X a 


M^i^ b 




Mmol 
(106 Mq) 










■ 2 ■ 

Sin I 
(106 Mq) 


[CO/H2] 

(106 Mo) 




Mo 




He 2-10 


11 


1.0 


330 


1.0 


140 


0.42 


475 


1.4 


NGC 5253 


15 


7.0 


94 


0.004 


37 


0.39 


680 


2.5 


NGC 1569 


50 


5.7 


26 


0.007 


4.6 


0.18 


140 


3.1 


NGC 3077 


50 


0.9 


93 


0.04 


4.4 


0.05 


45 


8.6 


Haro 2 


6 


3.5 


1100 


2.1 


750 


0.68 


2300 


1.1 


Haro 3 


50 


3.5 


430 


0.25 


110 


0.26 


790 


1.7 


H Zw 40 




6.1 






< 0.4 




<830 




Mrk 86 




7.7 






< 0.3 




<340 





"^Using the single-dish X vs. metallicity relation of Arimoto et al. 1996: log {^ x^ i ) ~ ~[^/-^] + ^-^"^ 
^Using: M^i^ = 210 Rpc AV^^^ Mq, (MacLaren, Richardson & Wolfendale 1988) 
''Mga, = Mmoi + Mhi, (see Table 1) 
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